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Abstract 

We have investigated analytically quantum tunneling of large spin 
in the biaxial spin system with the magnetic field applied along the 
hard and medium anisotropy axes by using a purely quantum-mechanical 
approach. When the magnetic field parallels the hard axis, the tunnel 
splittings of all the energy level pairs are oscillatory as a function of 
the magnetic field. The quenching points are completely determined 
by the coexistence of solutions of the Ince's equation. When the mag- 
netic field points the medium axis, the tunnel splitting oscillations 
disappear due to no coexistence of solutions. These results coincide 
the recent experimental observations in the nanomagnet Fes- 

In recent years much attention has been paid to quantum tunneling of 
magnetization (large spin) in nanomagnets, both from experiment and from 
theory [1]. Several magnetic particles have been identified as promising can- 
didates for the observation of such macroscopic quantum phenomena, where 
the magnetization (or the Neel vector) tunnels from one potential minimum 
to another one. The excellent examples that have widely studied are the 
molecular nanomagnets Fe§ [2-5] and Mn 12 [6-8] , which have the well-defined 
structures and magnetic properties. On the one hand, these phenomena are 
very interesting from a fundamental point of view because they extend our 
understanding of the transition between quantum and classical behavior. On 
the other hand, tunneling of the magnetization changes the magnetic prop- 
erties of the nanomagnets, which has the potential application for the data 
storage technology, e.g., making qubits - the elements of quantum computers. 

Recently Wernsdorfer and Sessoli [9] have observed a novel phenomenon, 
i.e. oscillations of tunnel splittings of the ground and excited states in the 
nanomagnet Feg described well by the spin Hamiltonian [2-5] 

H = AS 2 X - BS 2 y - gfi B S ■ H, (1) 
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where S is a spin operator, H is the magnetic field applied in the x-z plane, 
the spin quantum number s = 10, A pa 0.092K, B pa 0.229K, g pa 2 is the 
(/-factor, n B is the Bohr magneton. The zero-field Hamiltonian has a biaxial 
symmetry with hard, easy and medium axes along x, y and z, respectively. 
When H rotates from x to z direction, the oscillations of tunnel splittings of 
all the level pairs gradually disappear. 

In fact, oscillation of the ground state tunnel splitting AE S of the model 
(1) with the magnetic field H x (H || x) was predicted by using instanton 
technique [10]. The tunneling of spin is quenched when 

n 1 

h = h*(l -), n = 0,l,...,2s-l. (2) 

Here, h = % L , H c = 2s ( A+B ) j s the critical field at which the energy barrier 
h c ' c g^B toJ 

vanishes, —h*<h<h* = \f-^r^, and s is an integer or half odd integer. 
This kind of topological quenching is the result of quantum interference of 
different instanton paths within the context of macroscopic quantum tunnel- 
ing [11, 12], and need not be related to Kramers's degeneracy. Up to the first 
order of the rate -^g, the formula (2) was rederived by quantum- mechanical 
perturbation theory [13]. Very recently, Garg extended his previous work [10] 
to the excited states by using a discrete Wentzel-Kramers-Brillouin approach 
[14]. To order s" 1 , the quenching points for the excited state pairs are the 
same as those for the ground state pair. With increasing the magnetic field, 
the quenching points gradually decrease and finally disappear. These results 
were also obtained by the potential field description of spin systems with 
exact spin-coordinate correspondence [15] and the instant technique [16, 17]. 
Because quantum tunneling of spin in the nanomagnets was observed at very 
low temperatures, the quantization of spin levels becomes very important to 
explain well the experiments. In essence, the energy spectrum of the spin 
systems can help us to understand the mechanism of spin tunneling thor- 
oughly. In this paper, we have analytically diagonalized the Hamiltonian (1) 
in the framework of the Schrodinger's picture of quantum mechanics. The 
energy spectrum of the spin model is obtained in the large-s limit. It is 
clearly shown that when H || x, the tunnel splittings of the ground and ex- 
cited state pairs are oscillatory as a function of H x and the quenching points 
agree with the numerical simulation of model (1). When H || z, the tunnel 
splitting oscillations of all the energy level pairs disappear, which coincided 
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with those obtained by the phase space path integral [18]. These phenomena 
have also been observed in the experiments [9, 19]. 

Let E and $ m be the eigenenergies and eigenstates of H, respectively, 
then the eigenvalue equation in the basis \s,m > z reads 

«m-l$m-2 + U m+1 ® m+2 - t m _i $ m -l - t m+ l$ m+1 
1 



where 



-{-E + ^{A - B)[s{s + 1) - m 2 ] - gix B H z m}<$> m = 0, (3) 
u m ±i = \(A + + l)-(m±l)2]2-(m± 1)2, 



1 / 11 

*m±| = 2 9 ^ bHx ^ s ( s + !) - ( m - + 4 

. Obviously, it is very difficult to strictly solve the equation (3) for arbitrary 
s. However, in the large-s limit, Eq. (3) becomes [20] 

2^ 2$ o d ® r ^ 111 



cfe 2 da; L A + fi 4 4 1 - x 
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As(s + 1)^ 2N gn B H z ^Js(s + l) 

H i ^— (1 — # ) — ^ x 

A + B v ; A + B 



gfj, B H x y/s(s + 1) 



A + 5 
- th 

deriving Eq. (4), we have used 



Vl-x 2 ]$ = 0. (4) 



Here, x = , m and only the leading terms are remained, i.e. 0(s ). In 

y/s(s+l) 



<$>' S>" 
rF> „ lrf\ _l_ mil i mil , 

0(s + 1) 2s(s + l) 

1 d , _ 1 d 2 , ^ x 

0(s + 1) 2s(s + l)dx 2 

1 r d& , 1 d , d§. . 

0(s + l) dx 0( s + i)dx dx 
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u d 2 $ 

+ 



2s(s + 1) dx 2 

<]>' 



where 



u 



2 2y/ S (s+l) 
= t*± . 1 j-(t$) + --- 

2^/s(s + l)^ 
2^/s(s + l) dx 



\{A + 5)[ v /s 2 (s + 1) 2 (1 - x 2 ) 2 - s(s + l)x 2 
\{A + B)[s(s + 1)(1 - x 2 ) - x 2 /[2{l - x 2 )} 



and t ~ \g^BH x ^l s(s + l)Vl — ^ 2 for large s. 

Taking the transformations: $ = (1 — x 2 )~^y{x) and x = sin(2t), and 
substituting them into Eq. (4), we finally obtain the Hill's equation [21] 



d 2 v b 2 

— f + [A + (a - 6)cos(2t) + (c - 6)sin(2i) + — cos(At)]y = 0, (5) 
dt A o 



where 



-AE + 2As(s + 1) , Agfi B H x Js(s + l) 

A = : ^ , a = o v . 

A + B ' A+fi 



fc = ±4 + D + 

V a + 5 ' A + fi v ; 

Up to now, we have mapped the spin problem (1) onto a particle problem 
(5). The energy spectrum of Ti is completely determined by the characteristic 
levels of the Hill's equation. We note that when A = and H = 0, A = m 2 . 
So E = —B(y) 2 are nothing but the eigenvalues of the Hamiltonian (1) with 
integer s for even m or with half odd integer s for odd m. 

For Eq. (5), there exist two monotonically increasing sequences of real 
numbers a , an-, b 2 i, o! 2i _ x and b' 2i _ x {i = 1,2, •••) such that Eq. (5) has a 
solution with period n if and only if A = a , a 2 « or 6 2i , and a solution with 
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period 2n if and only if A = a' 2i _ x or b' 2i _ x . The ao,a2i,b2i,a' 2i __ 1 and b' 2i _ x 
satisfy inequalities: 

a o < b[ < a[ < b 2 < a 2 < b' 3 < a' 3 < 6 4 < a 4 < • • • . (7) 

According to the relation between E and A in Eq. (6), it is easy to see that 
the ground state of 7i corresponds to the allowed highest characteristic level 
of the Hill's equation with period 7r or 2n, depending on integer or half odd 
integer spin s. The lower the characteristic level of Eq. (5) is, the higher the 
associated eigenstate of the Hamiltonian (1) is. The tunnel splitting of the 
characteristic level pair of Eq. (5) with period 7r is 

AA 2i = a 2i - b 2i , (8) 

and the tunnel splitting of the characteristic level pair of Eq. (5) with period 
2ir is 

^2i-l — a 2i-l ~ Ki-i- (9) 
The tunnel splitting of the energy level pair of the Hamiltonian (1) can be 
evaluated by the tunnel splitting of its associated characteristic level pair 
using Eq.(6). We note that when two of three parameters A, H x and H z 
vanish, the Hill's equation (5) becomes the well-known Mathieu equation 
[22] . So the tunnel splittings of all the energy level pairs of the Hamiltonian 
(1) are monotonous increasing rather than oscillatory in the three cases [20]. 
For arbitrary parameters A, B, H x and H z , it is difficult to solve analytically 
Eq. (5). Here we only consider the following two special cases. 

(i) H z — (i.e. c = b). In this case, Eq. (5) reduces to the Ince's 
equation, which has been studied in a set of literature due to its physically 
basic importance [21]. To observe the oscillations of tunnel splittings of 
model (1), we must find some vanishing points at which its eigenstates are 
degenerate, i.e. a 2i = b 2i or a^i-i = ^2?-i- This is equivalent to find the 
coexistence of solutions of the Ince's equation, which means that there exist 
two linearly independent solutions (one even and one odd) with period 7r or 
2tt. Very fortunately, due to the positive coefficient of the last term cos(4t), 
the Ince's equation has the coexistence of solutions [21] with period n when 

a = -2nb (10) 

and with period 2n when 

a = -(2n-l)b, (11) 
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where n = 0, ±1, ±2, • • •. Obviously, the quenching points of tunnel splittings 
of the characteristic level pairs with period n are exactly shifted by half a 
period relative to those with period 2ir. Under the conditions (10) and (11), 
the Ince's equation becomes the Whittaker equation [21] 

1| + [ A _ p bcos(2t) + jCos(4t)]y = 0. (12) 

Here, p = 2n + 1 or 2n for Eq. (10) or (11), respectively. For Eq. (12), when 
6^0, then a — > 0, a 2i and b 2i — > (2i) 2 , and a! 2 i-i and b' 2i _i — > (2i — l) 2 [23]. 

Integer spin s. Due to Eq. (6), the eigenstates of the Hamiltonian (1) 
correspond to the characteristic levels of the Ince's equation with period ir, 
i.e. E« = -\{A + B)a + E , Ef = -\{A + B)a 2i + E and E\ = -\{A + 
B)b 2 i + E , E = ^As(s + 1), i — 1, 2, • • • , s. So the tunnel splitting of the 
energy level pair (Ef, Ef) is 

AEi = E\ - Ef = ±(A + B)AA 2l . (13) 
From Eq. (10), we have AEi = when 



H. = (2 " + 1) V- 4(A + B> . (14) 
Mb 

The existence of the quenching fields (14) clearly shows that the tunnel split- 
tings of all the energy levels of Ti are oscillatory as a function of H x and the 

period of oscillations AH = 2 V A ( A+B 1^ T^jg coincides with Garg's result (2) 

1 9/J-B to v ' 

[10, 14]. For Eq. (12), when \p\ = 21 + 1, then the even intervals of instability 
on the A axis disappear with at most I + 1 exceptions [21]. In other words, 
the characteristic values ao, a 2 i and b 2 i satisfy 

a < b 2 < a 2 < ■ ■ ■ < b 2i < a 2 i 

< b 2 (i + i) = a 2 (i+i) < ■ ■■ <b 2s = a 2s . (15) 

This means that the tunnel splittings of the / highest excited state pairs of the 
Hamiltonian (1) do not vanish while those of the other s — l energy level pairs 
vanish. With increasing the magnetic field H x (i.e. the quenching points 
gradually decrease and finally disappear when I > s. The configuration of 
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the quenching points agree with that from the numerical simulation of the 
Hamiltonian (1) (see Fig. 1). 

Half odd integer spin s. The eigenstates of the Hamiltonian (1) correspond 
to the characteristic levels of Eq. (12) with period 2ir, i.e. E°-\ = — \{A + 

B)a' 2i _ 1 + E and^'i = -\{A+ B)^^ + E ,i = l,2,---,s + ±. The tunnel 

splitting of the energy level pair (E a \ , E b '_ 1 ) reads 

1 2 1 2 

AE'^ = E*_i - El h = \{A + B)AA' 2i _ v (16) 
According to Eq. (11), we obtain A£'_i = when 

, 2nJA(A + B) 
K = — ^ (17) 

9Vb 

Obviously, the period of oscillations of tunnel splittings for all the energy 

level pairs (E a ' i,E b ' i) is AH' = W A ( A+B \ ^jch is the same with that of 
1 v 1-5' %-\> gvB ' 

the energy level pairs (E^,E b ). However, the quenching points for half odd 

integer spin s are shifted half a period (i. e. ^ A g^ B -) relative to those for 
integer spin s. For Eq. (12), if \p\ = 21, then at most I + 1 odd interval of 
instability on the A axis remain, i. e. 

b[ < a[ < b' 3 < 4 < • • • < U 2l _ x < a'a.! 

< K1+1 = a'a+i < • • • < b 2s = a f 2s- (18) 

It is easy to see that the tunnel splittings of the s — I + \ lowest energy level 
pairs vanish, but the tunnel splittings of the other / highest energy level pairs 
do not. When I > s — |, there do not exist the quenching points. These 
results also coincide with those obtained by the discrete WKB approach [14]. 

(ii) H x = (i.e. a = b). Let t — > t + |, then Eq. (5) becomes the Ince's 
equation 

^| + [A + ( c - 6)cos(2t) - ^cos(4t)]y = 0. (19) 

Because the coefficient of the last term cos(4t) is negative, Eq. (19) does not 
possess the coexistence of solutions, i.e. two linearly independent solutions 
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[21]. This means that the sign of equality in the inequalities (7) cannot hold 
and each energy level of 7i is singlet in the parameter space. So the tunnel 
splittings of all the level pairs are not oscillatory with H z , which coincide the 
experiment [9, 19]. 

In conclusion, we studied quantum tunneling of large spin in the biaxial 
spin systems by using quantum mechanics. The energy spectrum of the spin 
model (1) is obtained by solving Hill's equation (5), which is derived directly 
from the eigenvalue equation of the spin problem in the large- s limit. It is 
surprising that when H || x, the vanishing points of the tunnel splittings of 
all the energy level pairs obtained here by the coexistence of solutions of the 
Ince's equation coincide those given by the WKB method [10, 14] and other 
approaches [15-17]. However, our approach is a natural way of explaining 
the oscillations of tunnel splitting in the biaxial spin system, which is also 
applied in other spin Hamiltonians. 
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Fig. 1 The exact quenching points of all the level pairs of the Hamiltonian 
(1) with s = 10,A = 0.092/T, B = 0.229^, and H z = 0. 
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